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ABSTRACT 

Study of the primary anisotropics of the Cosmic Microwave Background (CMB) can be used to deter- 
mine the cosmological parameters to a very high precision. The power spectrum of the secondary CMB 
anisotropics due to the thermal Sunyaev-Zel'dovich Effect (SZE) by clusters of galaxies, can then be 
studied, to constrain more cluster specific properties (like gas mass). We show the SZE power spectrum 
from clusters to be a sensitive probe of any possible evolution (or constancy) of the gas mass fraction. 
The position of the peak of the SZE power spectrum is a strong discriminatory signature of different 
gas mass fraction evolution models. For example, for a flat universe, there can be a difference in the I 
values (of the peak) of as much as 3000 between a constant gas mass fraction model and an evolutionary 
one. Moreover, observational determination of power spectrum, from blank sky surveys, is devoid of any 
selection effects that can possibly affect targeted X-ray or radio studies of gas mass fractions in galaxy 
clusters. 

Subject headings: Cosmology: observations, cosmic microwave background; clusters : intracluster 
medium, gas fraction 



1. INTRODUCTION 

Clusters of galaxies, being perhaps the largest gravita- 
tionally bound structures in the universe, are expected to 
contain a significant amount of baryons of the universe. 
Moreover, due to their large angular sizes, observational 
estimates of their total mass My, the gas mass Mg and 
hence the gas mass fraction {fg = ^^) are easier. These 
estimates can be used as probes of large scale structure and 
underlying cosmological models. For example, the cluster 
fg would give a lower limit to the universal baryon fraction 
Jlfc/rim- Determination of fg has been done by numerous 
people (White & Fabian, 1995; Mohr et al. 1999; Sadat & 
Blanchard, 2001) and the values are in agreement within 
the observational scatter. A point to be noted here is that 
the estimated fg depends on the distance to the cluster 

(i.e fg oc dang )• Hcucc, if fg is assumed to be constant, 
then in principle, one can use the 'apparent' evolution of 
fg over a large redshift range to constrain cosmological 
models (Sasaki 1996 ). 

The question as to whether there is any evolution (or 
constancy) of gas mass fraction, however, is still debat- 
able, with claims made either way. For example, Schindler 
(1999) has investigated a sample of distant clusters with 
redshifts between 0.3 to 1 and conclude that there is no 
evolution of the gas mass fraction. Similar conclusion has 
been drawn by Grego et al. (2000). On the contrary, Ettori 
and Fabian (1999) have looked at 36 high-luminosity clus- 
ters, and find evolution in their gas mass fraction (in both 
sCDM and ACDM universes). See also David et al. , 1995; 
Tsuru et al. 1997; Allen & Fabian, 1998; Mohr et al. 1999). 
Observations suggest that, though fg of massive clusters 
{Tf,> 5KeV) appears to be constant, low mass clusters 
with shallower potential wells may have lost gas due to 
preheating and/or post-collapse energy input (David et 



al. 1990, 1995; Ponman et al. 1996; Bialek et al. 2000). It 
is also well known that ICM is not entirely primordial and 
there is probably continuous infall of gas, thereby, increas- 
ing fg with time. Thus, there is considerable debate re- 
garding the evolution of gas mass fraction. 

The intracluster medium (ICM) has been probed mainly 
through X-ray observations, but also through the so called 
Sunyaev-Zeldovich effect (Zel'dovich & Sunyaev, 1969) in 
the last decade (see Birkinshaw 1999 for a review). The 
SZE effect from clusters is a spectral distortion of the CMB 
photons due to inverse Compton scattering by the hot ICM 
electrons, with its magnitude proportional to the Compton 
y-parameter, given by ?; = -^f^ / UeTedl. Here, ks is the 
Boltzman constant, ctt is the Thomson scattering cross 
section, TOe is the electron mass, and rig and Tg are the 
ICM electron density and the temperature. Using bolo- 
metric SZE measurements, the fg has been obtained for a 
number of clusters (see Birkinshaw 1999). Recently, Grego 
et al. (2000) have made interferometric observations of SZE 
from a sample of 18 clusters. A major advantage of SZE 
over X-ray measurements is, SZE does not suffer from the 
(1 + z) * 'cosmological dimming', which makes SZE an 
useful probe of evolution of cluster gas mass fraction. 

Other than the targeted SZE observations, non-targeted 
'blank sky' surveys of SZE are one of the main aims of 
future satellite and ground based small angular scale ob- 
servations (Holder k Carlstrom, 1999; or with AMiBA ). 
Once the power spectrum is extracted from observations, 
comparison can be made with theory, to constrain cosmo- 
logical parameters and relevant cluster scale physics. The 
SZE power spectrum as a cosmological probe has been 
well studied (Refregier et al. , 1999; Komatsu & Kitayama, 
1999), although its use as a probe of ICM has seldom been 
looked at. 

Keeping such surveys in mind, in this Letter, wc look at 
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the SZE power spectrum as a probe of the ICM. We show 
it to be a very sensitive probe of the evohition (or con- 
stancy) of fg . Measurements of the primary anisotropy 
would give us 'precise' values of cosmological parameters 
(like ft,, r2„i, r^A, fi;,)- Hence, for our calculations, we as- 
sume that we know the values of cosmological parameters 
and do not worry about their effect on the SZE power 
spectrum. Any feature of the SZE power spectrum is at- 
tributed to specific cluster physics (like gas content). We 
note that this method of probing fg is not biased from any 
selection effect that can occur while doing pointed SZE ob- 
servations of X-ray selected clusters of galaxies, and hence 
is more desirable. 

Current observations of primary CMB anisotropies sug- 
gest a flat universe with a cosmological constant (Padman- 
abhan and Sethi, 2000). For our calculations, we take a 
flat universe with 0.^ = 0.35, fif, = 0.05 and h = 0.65 as 
our fiducial model. 

The paper is structured as follows. In the next section 
we discuss the distribution of clusters and model the clus- 
ter parameters. In §3, we compute the Poisson and clus- 
tering power spectrum from SZE and, finally, we discuss 
our results and conclude in S4.. 



2. THE SUNYAEV-ZEL'DOVICH EFFECT POWER 

SPECTRUM 

2.1. Distributing the galaxy clusters 

We set up an ensemble of galaxy clusters with masses 
between 10^^ < M < 10^^ M0, using the abundance of 
collapsed objects as predicted by a modified version of 
the Press-Schechter (PS) mass function (Press & Schechter 
1974) given by Sheth and Tormen (1999). The lower mass 
cutoff signifies the mass for which one expects a well de- 
veloped ICM. The gas mass is supposed to sit in the halo 
potential and arc distributed in the same manner. Note, 
the results arc not sensitive to either lower or upper cutoff. 
We probe up to rcdshifts of 5. Most of the power, however, 
comes from objects distributed at z<l.. 

We use the transfer function of Bardeen et al. (1986), 
with the shape parameter given by Sugiyama (1995) and 
the Harrison - Zel'dovich primordial spectrum to calcu- 
late the matter power spectrum Pm(k). The resulting 
COBE - FIRAS normalised (Bunn & White, 1997) mass 
variance (erg) is 0.9 for our fiducial model. 



2.2. Modelling the cluster gas 

We assume the ICM to follow a /3-profile with /3 = 2/3 

for simplicity. The other physical parameters of the clus- 
ters are determined using the virial theorem and spher- 
ic;al collapse model. We closely follow Colafrancesco & 
Vittorio (1994) in our modelling. We have for the gas 

/ 2 X -3/3/2 

density, ne(r) = Ue.o ( 1 + 72 ) • We take the gas to be 
extended up to i?„ = pr^ with p = 10. The central gas den- 
sity, ne,o is given by ne,o = fg m^{i°+x) where X = 0.76 is 
the average proton mass fraction and po is the central gas 
mass density. To account for the fact that there is a final 
cutoff in the gas distribution we introduce a Gaussian fil- 
ter at the cluster edge Ry given by ne(r) — > ne{r)e~^ /i^v^ 
where ^ = 4/7r is the fudge factor. 



We parametrize the gas mass fraction as 



fg = fg0{l + Zy 



M 



(1) 



where the normalization is taken to be fgo = 0.15, is based 
on local rich clusters. Wc look at combinations of both 
mass and redshift dependence for a range of evolutionary 
models. In particular, we look at a case of strong evo- 
lution given by fc = 0.5, s = 1 (Colafrancesco & Vittorio, 
1994); k — 0.1, s — 0.5 (as suggested by Ettori and Fabian, 
1999); k = 0.1, s = 0.1 (weak evolution); fc = (no mass 
dependence) and s = (no redshift dependence). 
For the core radius Tc and the temperature, we use 



,^ ,^ , imh-^Mpc 1 



1 + z 



(2) 
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fcsTe = 7.76/J- 



M 



2/3 



(1 + z)KeV. (3) 



Here, Ac(z) is the cluster overdensity relative to the back- 
ground. 

Putting everything in, we have the temperature distor- 
tion to be = g{x)y{9), with 



rUeC^ 



(4) 

The angular core radius 9c = rc/dang- The spectral form 
of the thermal SZ effect is given by 



9{x) = 



x^e^ 



(e^ - 1)' 



oth{x/2) - 4] , 



(5) 



where x = hv/ksTcmb- This specific spectral dependence 
of the thermal SZ effect can be used to separate it out 
from other CMB anisotropies (Cooray et al. , 2000). 

3. COMPUTING THE POWER SPECTRUM 

The fluctuations of the CMB temperature produced 

by SZE can be quantified by their spherical harmonic 
coefficients aim, which can be defined as AT(n) = 
Tq~^ X^im 0-imYimi^)- The angular power spectrum of SZE 
is then given by C; =< |a;m|^ >, the brackets denoting 
an ensemble average. The power spectrum for the Pois- 
son distribution of objects, can then be written as (Cole 
& Kaiser 1988, Peebles 1980) 



Poisson 



Jo dZ JMrni 



dM ' 



(6) 

where V{z) is the comoving volume and dn/dM is the 
number density of objects. 

Since these fluctuations occur at very small angular 
scales, we can use the small angle approximation of Leg- 
endre transformation and write yi as the angular Fourier 
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transform of y{0) as yi = 2n J y{0)Jo[{l + l/2])0]ede (Pee- 
bles 1980, Molnar & Birkinshaw 2000). 

In addition to Poisson power spectra, one would expect 
contribution to a 'correlation power spectrum' from the 
clustering of the galaxy clusters. Following Komatsu and 
Kitayama (1999), we estimate the clustering angular power 
spectrum as 

^Clustering ^ r"""^ rf^^Kifl X (7) 

Jo dz 

r n 2 

r^^' rfM ^^^^f"^ (M, z,MM, z) , 

where b{M, z) is the time dependent linear bias fac- 
tor. The matter power spectrum, Pm{k,z), is related 
to the power spectrum of cluster correlation function 
Pc(fc, Ml, M2, z) through the bias, i.e Pc{k, Ml, M2, z) = 
b{Ml,z)b{M2,z)D'^(z)Prn{k,z = 0) where we adopt 
b{M,z) given by 6(M, z) = (1 + 0.5/i^'*)"-"^-"-"2"(l + (i/2 _ 
l)/(5c) (Jing 1999 for details). This expression for the bias 
factor matches accurately the results of N-body simula- 
tions for a wide range in mass. In the above equation D{z) 
is the linear growth factor of density fluctuation, 6c = 1-68 
and u = Sc/cr{M). 

4. RESULTS AND DISCUSSIONS 

Wc study the power spectrum of SZE from clusters of 
galaxies, under the assumption of a 'precise' and 'a pri- 
ori' knowledge of the cosmological parameters. We also 
assume that in the Z-range of relevance, thermal SZE from 
clusters of galaxies are the dominant contributors to the 
temperature anisotropy. The other secondary anisotropics 
are either smaller in strength or contribute at even higher 
I's or have different spectral dependence (Aghanim et al. , 
2000; Majumdar et al. , 2000). 

We have plotted the Poisson SZE power spectrum in 
Figl, left panel. Clearly, the primary feature distinguish- 
ing a non-evolutionary constant fg model from an evolu- 
tionary one is the position of the peak. The model with a 
constant fg peaks at a higher Z-value and also has greater 
power. The constant fg model peaks at / ^ 4000. This 
result is in agreement with that of Komatsu & Kitayama 
(1999). If one assumes that there is no evolution of fg 
with redshift (i.e s=0), the peak is at Z ~ 1100, whereas in 
the case of no dependence on mass (k=0), the peak is at 
/ ~ 2500. Based on EMSS data (David et al. , 1990), Co- 
lafrancesco & Vittorio (1994) (and also Molnar & Birkin- 
shaw 2000) model fg with k=0.5 & s=l. For this case, 
we see that the turnover is at a very low I ~ 900. As- 
suming a mild evolution (k=0.L s=0.1), we get the peak 
at I ~ 2100. We also show results for (k=0.5, s=0.5) and 
(k=0.1, s=0.5). The last parametrization is based on the 
recent analysis of ROSAT data by Ettori & Fabian (1999). 
It is evident that the difference, in the l-vahie of the peak 
of the constant fg scenario from an evolutionary one, can 
range between I ~ 1500 — 3200.. The position of the peak 
thus is a strong distcriminatory signature of any evolution 
of/,. 

It is easy to understand the shift in the peak of the SZE 
power spectrum. Let us consider the case s=0, i.e. fg 
depends only on total mass. From Eqn.(l), this means 
an enhanced reduction of fg of smaller mass clusters rel- 



ative to the larger masses and so a reduction of power 
at larger I's, (since smaller masses contribute at larger I). 
Hence, the peak shifts to a lower I. For the case k=0, 
(i.e only redshift dependence), we now have structures at 
high z contributing less to the power (than without a red- 
shift dependence). Since from PS formalism, less massive 
structures are more abundant at high z, this negative de- 
pendence of fg on redshift cuts off their contribution more 
than the more massive clusters. Hence, once again there 
is less power at high I and the peak shifts to lower Z-value. 
The parametrization of Eqn. (1) affects the larger masses 
less, as evident from almost equal power seen at l< 600, 
for all models. The net effect is a reduction of power at 
smaller angular scales, and hence a shift in the position of 
the peak to a smaller multipole value. 

We note that, these results are irrespective of the argu- 
ments given (see Rines et al. , 1999) to explain any possible 
evolution of fg . In their case, they assume fg to be con- 
stant and relies on the cosmology to change the angular 
diameter distance, so that there is an 'apparent' change 
in fg . In such a case, if there is 'actually' even a slight 
evolution of fg , then one can still account for it with a 
non-evolutionary model, by simply changing the cosmo- 
logical parameters. Our method does not assume a priori 
any constancy (or evolution) of fg and tries to look for it. 

In Figl, right panel, we show the SZE clustering power 
spectrum. For all models, it falls of at a smaller I w.r.t 
Poisson power spectrum. Since for clustering, the peak 
depends on the average inter-cluster separation, which is 
fixed once the cosmology is fixed, there is no appreciable 
spread of the peaks in Z-space. The only difference is in 
their relative power w.r.t each other which depends on the 
total gas mass available to distort the CMB. Addition of 
the clustering power spectrum to the Poisson case results 
in slight shift of the peaks to lower Z's. 

It maybe possible to measure the power spectrum of 
SZE with the ongoing and future high angular resolution 
CMB observations. In principle, observations with SUZIE, 
OVRO, BIMA and ATCA can probe the range in I from 
» 1000 - 7000 and a frequency range of « 2 - 350 GHz. 
The SZE power spectrum would also be measured with 
increased precision by the proposed ALMA and AMiBA 
(which is geared for blank sky surveys). 

Finally, let us comment on the validity and robustness 
of our results. In Fig2, we show results for an open uni- 
verse {Hq — 0.35, h = 0.65). It is clearly seen that the 
difference in the peak position of constant fg and evolu- 
tionary models remain far apart (Infact, for same param- 
eters of A; = 0.5, s = 1, the difference increases to » 4500 
from that of « 3000 in a flat universe). It is seen that 
the turnover of the SZE power spectrum is insensitive to 
the mass cutoff, since main contribution to the anisotropy 
comes from clusters with lO^'^M© < M < IO^^Mq. In 
Fig2, we also indicate the effect of having a more compact 
gas distribution with p= 7. We see that shift in the peaks 
are negligible (though the height is reduced a little) and an 
uncertainty as to how far the gas extends is not major. The 
use of a single /3 to model the full gas distribution intro- 
duces little error, though a /3-model fits the inner cluster 
regions better. This is because the major contribution to 
the anisotropy comes from around the core region, and in- 
creasing /3 slightly decreases the overall distortion, without 
touching the peak. Also, a modified M — T relation (more 
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suitable for ACDM) does not change the conclusions of this 
paper (although amplitude of distortion slightly changes). 
For a more detailed analysis, however, one should take bet- 
ter observationally supported gas density and temperature 
profiles (see Yoshikawa & Suto, f999). These points will 
be discussed in greater detail in a future publication. 

In conclusion, we have computed the angular power 
spectrum of SZE from clusters of galaxies. We have shown 
the position of the peak of the power spectrum to bear 
a strong discriminatory signature of different fg models. 
One of the goals of arc minute scale observations of the 
CMB anisotropy is to measure the SZE power spectrum 
from blank sky surveys. Such observational results can be 
used to constrain fg models. This also has the added ad- 
vantage of being devoid of uncertainties that can creep 



in through 'selection biases' in estimating the fg using 
pointed studies of X-ray selected galaxy cluster. Our 
method, thus, provides a powerful probe of evolution (or 
constancy) of gas mass fraction and can potentially resolve 
the decade long debate. 
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tacharya for valuable suggestions and Pijush Bhattachar- 
jee for encouragements. He also wishes to thank the Ra- 
man Reseach Institute, Bangalore. After the completion 
of this work, we learnt about a complimentary work by 
Carlstrom & Holder (2001) which looks at changes in gas 
profiles due to preheating and its effect on SZE. 
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Fig. 1. The Poisson (left panel) and elustering (right panel) power spectra due to SZE from galaxy clusters for different fg models. 
For both the panels, the thick solid line correponds to constant fg model, the thick dashed line has no evolution with redshift and the thick 
dash-dotted line has no evolution with total mass. The thin lines are for the cases: a) k=0.5, s=l; b) k=0.5, s=0.5; c) k=0.1, s=0.5 and d) 
k=0.1, s=0.1. 
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Fig. 2. — The Poisson SZE power spectra are plotted for different eosiiiologies and with different extension of the gas mass. The soUd hnes 
are for a ACDM, with = 0.35, Qa = 0.65, h = 0.65, and the thin Unes are for OCDM with = 0.35, h = 0.65. The OCDM hnes have 
been multiplied by a factor of 10 in the plot. The solid and dashed lines are for gas mass extending upto lOrc, whereas the dash-dotted and 
the dotted lines are for extension upto 7rc. 



